Traveling time and traveling length in critical percolation clusters 
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We study traveling time and traveling length for tracer dis- 
persion in two-dimensional bond percolation, modeling flow 
by tracer particles driven by a pressure difference between two 
points separated by Euclidean distance r. We find that the 



minimal traveling time tmin scales as t„ 



which is 



and 

1.21 



different from the scaling of the most probable traveling time, 
i ~ r^'''^. We also calculate the length of the path correspond- 
ing to the minimal traveling time and find £min ~ ' 
that the most probable traveling length scales as £ ~ r 
We present the relevant distribution functions and scaling re 
lations. 
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The study of flow in porous media has many appli- 
cations, such as hydrocarbon recovery and ground-water 
pollution ||-||. Here we study an incompressible flow on 
two-dimensional bond percolation clusters |^ at critical- 
ity where fluid is injected at point A and recovered at 
point B separated from point A by Euclidean distance r. 
At time i = we add a passive tracer Q at the injection 
point 1^. We investigate the scaling properties of the 
distributions of traveling time, traveling length, minimal 
traveling time, and the length of the path corresponding to 
the minimal traveling time of the tracer particles. We find 
new dynamical scaling exponents associated with these 
distributions. 

Our first step is to calculate the pressure difference 
across each bond by solving Kirchhoff's law, which is 
equivalent to solving the Laplace equation. The velocity 
across a given bond is proportional to the pressure differ- 
ence across the bond; we normalize the velocities assum- 
ing the total flow between A and B is fixed, independent 
of the distance between A and B and the realization of 
the porous media 

We simulate the flow of tracers using a particle- 
launching algorithm (PLA) where a tracer particle 
starting from the injection point A travels through the 
medium along a path connected to recovery point B [pdl . 
The probability pij that a tracer particle at node i selects 
an outgoing bond (ij) is proportional to the velocity of 
flow on that bond; pij — Vij / '^j^ Vik, where the k sum- 
mation should be taken over all outgoing bonds, i.e., for 
Vik > 0. In this process, the time taken to pass through 



the bond (ij) is inversely proportional to the velocity of 
that bond, i.e., tij = 1/vij. 

We measure the distributions, P{i) and P{i), of the 
traveling time t and the traveling length £ between A and 
B for 10000 tracer particles for each realization. We sam- 
ple over 10000 different realizations with the two points 
A and B fixed. For each realization, we also find the 
minimal traveling time and the path which corresponds 
to the minimal traveling time to obtain P{tmin) and 
P{£min)- We run the simulation for system size L x L 
where L — 1000 ^ r, and find a well-defined region 
where the distributions follow the scaling form 



(1) 



where x denotes imin, tmin, ^ or t. The normalization 
constant is given by A^ ^ (x*)^^ and we find the scaling 
functions to be of the form f{y) ~ exp(— aajj/""^^). The 
maximum of the probability is at x* . Simulation shows 
that X* has a power-law dependence on the distance r, 

a:*~r'*-. (2) 

The exponents (j)x and dx are related by (l)x = ^/{dx — 1) 
I p^ . The scaling function / decreases sharply when x 
is smaller than x* . The lower cutoff is due to the fact 
that the traveling distance cannot be smaller than the 
distance r. 

The path which takes minimal time is not always the 
shortest path. However we find that the distribution 
of tmin coincides with the distribution of the chemical 
lengths between points separated by distance r studied 
in detail in Ref. p| . 

In Figs. 1(a), 2(a), and 3(a), we show the log-log plots 
of distributions P{tmin), P{t), and P{i), respectively. 
For different distances r = 4, 8, 16, 32, 64, and 128, we 
determine the characteristic size a;* as the peak of the 
distribution. In Figs, lb, 2b, and 36, we plot x* ver- 
sus distance r in double logarithmic scale and linear fit- 
ting yields the exponents dx for each distribution. In 
Figs. 1(c), 2(c), and 3(c) we collapse the data by rescaling 
X by its characteristic size x* . All distributions are con- 
sistent with the scaling form of Eq. (Q) . The measured 
values of scaling exponents are summarized in Table I. 

As shown in Fig. |^(b), the most probable traveling 
length i* scales as £* ~ r'^i where = 1.21 ± 0.02. Note 



1 



that di is significantly different from the minimal path 
exponent dmin = 1-130 ± 0.002 [|l^, while it is within 
the error bars of the exponent for the optimal path in 
random energy landscapes, dopt — 1-2 ± 0.02 ||l^, and 
the shortest path in invasion percolation with trapping, 
dopt = 1.22 ±0.01 

In many transport problems, the characteristic time 
scales with the characteristic length with a power law, 
t* ^ {ty . Since t* scales as r''* and I* scales as r'^*, it 
is reasonable to assume that z — dt/dg. Combining this 
relation, the relation i ~ Eq. ([T]), and the identity 
P{(-min)dt.min = P {tmin)dtmin, wc obtain scahng rela- 
tions between exponents, 

- l)d,„„. = (gt,„.„ - l)dt„.„ (3) 

This scaling relation is well satisfied by the set of scaling 
exponents given in Table I. 

Because of flow conservation, the velocity at distance 
r' from point A should scale inversely proportional to 
the number of bonds at this distance, which scales as 
{r'Y'^^^ where ds is the fractal dimension of the trans- 
port backbone. Then the traveling time for a particle to 
travel the distance r is given by 

~nr)-[^f'-r^^. (4) 

Note that t*{r) is the most probable traveling time in 
our system, so we obtain the scaling relation di = ds- 
Thus, the most probable traveling time is characterized 
by the transport backbone dimension of the media. This 
result is consistent with the homogeneous case, where 
ds = 2. The most recently reported value for the fractal 
dimension of the backbone is dg = 1.6432 ± 0.0008 |l|] 
for d — 2, which is in agreement with our results (Table 

I)- . . 

The minimal traveling time is the sum of inverse ve- 
locities over the fastest path where as noted above the 
fastest path is statistically identical to the shortest path. 
While the velocity distribution has been studied exten- 
sively (e.g. it is known to be multifractal), because 
the velocities along the path are correlated, how the min- 
imum traveling time distribution is related to the local 
velocity distribution is an open challenge for further re- 
search. 
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X 


dx 


9x 


f ■ 


1.13 ±0.01 


2.14 ±0.05 




1.33 ±0.05 


1.90 ±0.05 


£ 


1.21 ±0.02 


2.00 ±0.05 


t 


1.64 ±0.02 


1.62 ±0.05 



Table I. Results for the exponents. Our d = 2 results 
for and gi^^^ are within error bars of dmin and g'^ 

in Ref. |jl^. For comparison, the theoretical values of dx 
and gx for d — Q are all 2. 
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FIG. 1. (a) Log-log plot of the minimal traveling time 
distribution P{tmin) for separations r — 4,8,16,32,64, and 
128 between injection and recovery points, (b) Log-log plot of 
the most probable minimal traveling time versus r. A linear fit 
yields dt„i„ = 1.33±0.05. (c) The data obtained by rescaling 
the minimal time with its characteristic time f^in ~ r^'^^. A 
fit of the power-law regime gives = 1-90 ± 0.05. 
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FIG. 2. (a) Log-log plot of traveling distance distribution 
P{I) for r = 4, 8, 16, 32, 64, and 128. (b) Log-log plot of the 

most probable traveling length versus r. A linear fit yields 
di = 1.21 ± 0.02 (c) The data obtained by rescaling the trav- 
eling length with its characteristic length £* ~ r^'^^. A fit of 
the power-law regime gives = 2.0 ± 0.05. 



FIG. 3. (a) Log-log plot of traveling time distribution 
P{i) for r = 4, 8, 16, 32, 64, 128. (b) Log-log plot of the most 

probable time versus r. A linear fit yields df — 1.64±0.02. (c) 
The data obtained by rescaling the time with its characteristic 
time i* ^ r^'''"^. A linear fit of the power- law regime gives 
g{ = 1.62 ± 0.05. 
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